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Temperature and velocity oscillations have been found in a rectangular Rayleigh-Benard convec- 
tion cell, in which one large scale convection roll exists. At Ra = 8.9 x 10^^ and Pr — 4, temperature 
oscillation can be observed in most part of the system and the oscillation period remains almost 
constant, = 74 ± 2 seconds. Velocity oscillation can only be found in its horizontal component, 
Vy (perpendicular to large scale circulation plane), near the cell sidewall, its oscillation period is 
also constant, = 65 ± 2 seconds, at these positions. Temperature and velocity oscillations have 
different Ra dependence, which are respectively indicated by Peclect number Per = CSSTia" *^ and 
Pcv = Q.2diRa^'^^ . In comparison to the case of cylindrical cell, we find that the velocity oscillation 
is affected by the system geometry. 



INTRODUCTION 

The coherent large-scale flow structure is one of the 
most important constituents in turbulent thermal con- 
vections, which occur widely in the oceans, the atmo- 
sphere, the earth's outer core and galaxies |lj. Its de- 
tailed flow dynamics has been mainly investigated in 
the classical Rayleigh-Benard convection system, where a 
confined enclosure is heated from below and cooled from 
above 0, S 13- The system has three control parame- 
ters, namely the Rayleigh number {Ra), Prandtl num- 
ber (Pr) and the aspect ratio (F). The Rayleigh num- 
ber is defined as Ra — agl? jvK, with g being the 
gravitational acceleration, and a, v, and k being, re- 
spectively, the thermal expansion coefficient, the kine- 
matic viscosity, and the thermal diffusivity of working 
fluid. The Prandtl number is defined as Pr=: v jn,, and 
the aspect ratio F is the ratio of the lateral dimension 
of the system to its height. In the system, a prominent 
self-sustained low frequency oscillation of the measured 
local velocity and/or temperature has been generally ob- 
served in numerous experiments and numerical simula- 
tions i, i, i, S, [t,, 8, 9, m tJI 111 [13, Q E [11, [13, 




19|,|20|,|21|,|2^|23|,|2J,|25|,l2al27|. However, the ori- 
gin of this characteristic frequency still remains elusive 
although much effort has been made to tackle this issue 
in the past decades. 

Explanation of low frequency oscillations can be traced 
to the "bubble model" ^] and later the "plume model" 
14j . which were developed based on the presence and 
persistence of thermal structures (e.i?. thermals, plumes) 
[3| . It is assumed that a thermal boundary layer thickens 
in time by the upward diffusion of the temperature field. 
When a critical thickness is attained, the boundary layer 
will erupt and quickly expel all the buoyant fluid through 
plumes and thermals away from the boundary. The de- 
pleted boundary layer thus grows once more by diffusion 
and the process repeats in a periodic manner. Such a 
phenomenological picture has been observed in high Pr 
fluid with the shadowgraph technique ^28^ . In later exper- 



iments, it has been observed that the unique frequency 
is related to the large-scale circulation (LSC) in a helium 
gas cylindrical cell with F = 1 [1, [3, [1, [111 , and its value 
is more or less the same as the circulation frequency of 
the large convective flow. A coupled-oscillator model was 
proposed to explain the observed phenomenon, in which 
the plume clusters from one plate are triggered by those 
from the opposite plate and the overall process will be 
self-sustained !?| . A prediction of the model is that the 
oscillation period is twice of the cell-crossing time of the 
clustered plumes. This has been verified by later direct 
velocity measurements 17, 18, 23] . AH the above men- 



tioned models suggest that the thermal plumes and the 
thermal boundary layers play a dominant role in the gen- 
eration of the oscillation. 

In thermal convection, it has been found that, for the 
most part, the LSC has a preferred orientation in F = 1 
cylindrical cells and it oscillates azimuthally around this 
orientation 12, [H [13, [13, [13]. Recently, several new dy- 
namic processes, e.g. reorientation, cessation and so on, 
have been observed in the cylindrical cells [13, [H, [13], 
but it has been argued that such processes have little 
relevance to the azimuthal oscillation [ll| . Based on the 
phenomenon of horizontal oscillations of LSC [13], a new 
dynamical model has been proposed 13, 27]. This model 
suggests that the horizontal oscillation of large scale fiow 
is the main source of the low frequency oscillation in the 
system. Instead of being caused by the periodic plume- 
ejections from the boundary layers, this horizontal oscil- 
lation is intrinsic to the bulk dynamics. The periodical 
ejections of plumes may then be viewed as passive objects 
of the process. 

The LSC dynamical picture can be simplified if the 
large scale reorientations or cessations are not consid- 
ered, the horizontal oscillation of LSC will be analogous 
to a classical harmonic oscillator wandering around its 
origin. In this case, one may wonder whether the system 
geometry plays a role in the horizontal oscillations of the 
LSC. In most cases experiments and simulations have 
been performed mainly in a cylindrical cell of aspect ra- 
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tio one in which axis-symmetry holds and only one large 




13, Ei, [3 H [17 

At low Ra, the 



scale convective roll exists [6|. L 

3, Ei, M, M, m M, m, m,M, \M 

linear stability analysis predicted that one asymmetric 
convection roll occupies the whole bulk regime when the 
cell provides an axis-symmetric boundary condition [s^] , 
and the convection roll could be aligned in any direction 
around the circle [ssj . But in real systems there are in- 
evitably imperfections, initial instabilities and other ex- 
ternal factors, which lead to the evolution of one stable 
branch, with all other branches being unstable 34|. In 
numerical simulations, similar picture could be expected 
for the case of high Ra convection [sH . 

In this paper, we report turbulent thermal convec- 
tion experiments in a rectangular cell, in which the axis- 
symmetry is absent. Local velocity and temperature 
measurements at various locations were made over vary- 
ing values of Ra and at a fixed Pr. We focus on the 
oscillatory behavior of measured local velocity and tem- 
perature fields. The remaining sections of the paper are 
organized as follows. In Sec. II we describe the rect- 
angular convection cell and the velocity and tempera- 
ture measurements. In Sec. Ill, we present and discuss 
the experimental results. There we first show the atuo- 
correlation functions and histograms of the velocity and 
temperature signals. We also discuss the spatial distri- 
bution of the velocity and temperature oscillations. We 
then present the Ra-dependence of the periods of these 
oscillations. A brief summary is given in Sec. IV. 



EXPERIMENT SETUP 

The convection cell has a square vertical cross section, 
its length, width, and height are 81 x 20 x 81 (cm), 
as shown in Fig. [T] (a). With this geometry, the axis- 
symmetry does not exist, and the large-scale flow is ex- 
pected to be largely confined in the plane with the aspect 
ratio r = 1. The working fluid was water. Detailed de- 
scription of the apparatus can be found in Ref . [s^ . Here 
the coordinates are defined as follows: the origin coin- 
cides with the cell center, its x axis points to the right, 
the z axis points upward, and the y axis points inward. 
We measured the velocity and temperature signals at po- 
sitions marked by the solid lines in Fig. [1] (b) . The veloc- 
ity was measured by Laser Doppler Anemometry (LDA) , 
the measurement volume was approximately 0.12 mm 
in diameter and 1.1 mm in length _16i. Along the line 
(x, 0, 0), two-dimensional velocity (vy, Vz) was measured, 
and velocity {vx, Vz) was measured along the line (0, 0, z). 
At each position ^ 10^ or at least 5.2 hrs statistically 
independent samples were acquired. The temperature 
measurement apparatus was similar to that described in 
Ref.[37|. The local temperature was measured by us- 
ing a thermistor about 300 /im in size. The thermistor 
served as one arm of an AC Wheatstone bridge driven si- 
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FIG. 1: (a) Schematic drawing of the convection cell and the 
coordinates for the experiment; (b) The x — z cross section 
plane of the cell at y = (shown as the dashed plane in (a)). 
Inside the plane, the large scale flow is shown as the gray cir- 
cle. The measured positions are along the marked solid lines, 
the positions where the oscillation can be found from tem- 
perature signal are indicated by the line-combs, those from 
velocity by the line-diamonds. 



nusoidally at IkHz in frequency and 0.5 V in amplitude. 
The output of the bridge was first fed to a lock-in ampli- 
fier and then digitized by a dynamic signal analyzer at a 
sampling rate of 128 Hz, the recording time was 4.5 hrs 
to ensure good data statistics. Ra and Pr was kept at 
8.9 X 10^^ and 4.0 respectively unless stated otherwise. 



RESULTS AND DISCUSSION 

Autocorrelation function and histogram 

To test whether a low frequency oscillation exists 
we analyze the auto-correlation function, C{t) =< 
SX{t)SX{t + t) > of these measured signals, here 6X 
is the fluctuation of X normalized by its RMS value, X 
can be Vx, Vy, Vz and T and <> denotes the time av- 
erage. Examples of the auto-correlations of Vy, Vz and 
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FIG. 2: (a) Top to bottom curves: auto-correlation functions 
of Vy (shifted up by 0.2), Vz and temperature (shifted down 
by 0.2) at (385 mm, 0,0); (b) same as (a) but at (240 mm, 
0,0); (c) top to bottom curves: auto-correlation functions of 
Vx (shifted up by 0.2), Vz and temperature (shifted down by 
0.2) at (0,0, -385 mm); (d) same as (c) but at (0,0, -240 
mm). 
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FIG. 3: Histograms of (a) Vy at (385 mm, 0, 0), (b) Vz at (385 
mm, 0,0), (c) Vy at (240 mm, 0,0) and (d) Vz at (240 mm, 
0,0). Their fitting curves are plotted in solid lines from (a) 
two-peak Gaussian function, (b) one-peak Gaussian function, 
(c) and (d) one-peak exponential power function. 



T at (385 mm, 0,0) near the sidewall are shown in Fig. [2] 
(a). An oscillation can be found in Vy, whose intensity, 
represented by the magnitude of the second order peak, 
is as strong as that observed in the vertical velocity mea- 
surement near the sidewall of cylindrical cells 1^, 2^, 21 1 . 
But the vertical velocity Vz here shows no oscillation from 
its auto-correlation function, which differs from that in 



the cylindrical cells 3, 13, 21 1. The oscillation of tem- 
perature signal is much weaker than that in the cylin- 
drical cell [2l|, but the second order peak can be well 
resolved from the background noise fluctuation. The os- 
cillation periods obtained from Vy and T are denoted as 
ty and tx respectively. When the measurement moves 
away from the sidewall to the region of the cell, as shown 
in Fig. [5] (b) where the location is (240 mm,0, 0), the ve- 
locity oscillation can not be observed from both Vy and 
Vz components, whereas temperature oscillation can still 
be detected. In the measurement along the line (0, 0, z), 
two examples at (0, 0, —385 mm) and (0, 0, —240 mm) 
are shown in Fig. [2] (c) and (d) respectively, in which only 
temperature oscillation can be observed, the velocity one 
cannot be detected in either Vx or component. 

We now examine the histograms of these velocity sig- 
nals, some examples are shown in Fig. [3) Along the line 
(a;,0,0), there are two peaks in the histogram of Vy at 
location (385 mm, 0,0) (Fig. [3] (a)), and it can be fit- 
ted with a two-peak Gaussian function. The distribu- 
tion of Vy has two peaks which is consistent with the 
strong horizontal velocity oscillation shown in its auto 
correlation function (Fig. (Ha)). The histogram of Vz 
(Fig. [3)^b)) exhibits typical Gaussian- like distribution, it 



is also consistent with the property of its auto-correlation 
function (Fig.jDJa)). At location (240 mm, 0, 0), both ve- 
locity components Vy and Vz have similar single-peaked 
histograms (Fig. ^c) and (d)), they can be fitted by a 

stretched exponential function (~ exp[ 



(^^]), here 



the fitted /3 is 1.6 ±0.1. We also examined the histograms 
for both the Vx and Vz velocity components along the 
line (0,0, z), and find single-peaked Gaussian-like dis- 
tributions for positions close to the bottom plate and 
stretched-exponential distribution for locations far away 
from the plate. Such results are consistent with the obser- 
vations that no oscillation in their auto-correlation func- 
tions (Fig.^b), (c) and (d)). 

It is a surprise that the oscillation in the velocity field 
cannot be observed in its vertical component Vz- Similar 
finding has also been reported in the PIV measurement 
of the present system [36|, where no velocity oscillation 
has been seen in the large scale circulation plane. Such 
an observation is different from that in the cylindrical 
cell and does not support the "plume model" explana- 
tions 0, S, [11 [13, \m M, If the thermal 
plumes are ejected from the conducting plates in a peri- 
odical manner [l3,[2^ 28 1, the vertical velocity Vz should 
contain an oscillatory component, as the thermal plumes 
are driven by the buoyancy force. The absence of this 
oscillation suggests that temperature and velocity oscil- 
lations may have different origins and be controlled by 
different mechanisms. 
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Spatial distribution 

One can see that both temperature and velocity os- 
ciUations can be observed in the present system (al- 
though not in all positions). The positions where the 
oscillation can be found from T and Vy are marked in 
Fig. [Tfb). The temperature oscillation can be detected 
in the range x G {80 mm, 395 mm} along the line (x, 0, 0) 
and z G { — 130 mm, —387 mm} along the line (0,0,2:). 
There is no temperature oscillation in the central region 
of the cell. We examine the data in the cylindrical cell of 
Ref. [3^ and find the same result. However in the cylin- 
drical cell, it has been found that the temperature oscil- 
lation near the sidewall is much stronger than that near 
the bottom plate [l9| . Here we cannot distinguish where 
the temperature oscillation is obviously stronger from the 
auto correlations, as shown in Fig. [21 . 

From Fig. \V[h) , we can see that velocity oscillation can 
only be detected in its component Vy near the sidewall 
{x £ {250 mm, 395 mm}). Correspondingly, in these po- 
sitions two peaks have been observed in the histogram of 
Vy. The values of the peak positions in the histograms 
are shown in Fig. H) The difference between the two peak 
values becomes smaller as the measuring position moves 
away from the sidewall to the central region, and the two 
values merge into one at x around 250 mm, where the 
oscillation becomes undetectable in the velocity field, as 
illustrated in Fig. [T] No oscillation can be found in other 
measured positions, including x £ {0, 250 mm} along the 
line (a;, 0,0) and the whole line (0,0, z). The result is in 
sharp contrast to what is observed in the cylindrical cell 



23l |. In the cylindrical cell, the oscillation could be 
found everywhere in the bulk regime from the horizon- 
tal velocity measurement, and the oscillating intensity is 
stronger in the center than that in other positions [23|. 
Recently, the oscillation has also been explored in a rect- 
angular cell, where the horizontal cross section is square 
and the aspect ratio F = 4 [26]. The authors found that 
the oscillation stems from the periodical switching be- 
tween two corotating convective rolls. Therefore, the os- 
cillations observed in the present system should not share 
the same dynamic mechanism as that case because only 
one large scale roll exists here ;36]. Unlike the recent 
convection experiment in a rectangular box 

Next, we discuss the position dependence of the oscil- 
lation periods. It and ty, obtained from temperature and 
velocity fields repectively. As shown in Fig.[5l we can see 
that within experimental uncertainty the temperature os- 
cillation period (It) may be regarded as a constant for 
all the positions measured, = 74 ± 2 seconds for both 
the vertical and horizontal scans. The fact that tem- 
perature oscillation can be observed in most part of the 
measured locations with the constant oscillation period 
implies that it is a global characteristic of the convective 
flow and the oscillation period is a typical time scale of 



in 

E 
E 




200 300 
X (mm) 

FIG. 4: The peak values of Vy as function of x when the 
measured position moves from the cell center to the sidewall. 
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FIG. 5: The measured oscillation periods from the velocity 
Vy along the line {x, 0, 0) (solid squares) and temperature sig- 
nals along the lines (x, 0, 0) (open circles) and (0, 0, z) (solid 
circles) . 



the system. Oscillation of the velocity field can only be 
observed in a rather narrow region close to the sidewall, 
inside this region the turbulent intensity was found to 
be the strongest in the system fs^. The velocity oscilla- 
tion period (ty) also remains more-or-less constant for all 
the measuring positions, ty = 65±2 seconds. Note that 
ty is about 14% less than <t, such a difference between 
the two oscillation periods is beyond the experimental 
uncertainties for the respective measurements. 



Ra dependence 

Keeping constant Pr= 4.0, we also measured the os- 
cillation periods [tr and ty) of temperature and velocity 
fields over the range of Ra from 2 x 10^° to 8.9 x 10^^. 
In order to compare to the results in previous literatures 
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FIG. 6: The Ra-dependence of Peclet number. Pcv from the 
velocity oscillation (open circles) and its fitting 0.29Ra''^'^ is 
plotted as dash line; Per from the temperature oscillation 
(solid squares) and its fitting is plotted as solid 

line; Peclet number obtained from the rotational rate in the 
PIV measurement [36l ] is plotted here as dotted line. 
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of the Peclet number {PeT{v) = ^ 



is normalized in the form 



tT(,,) 



subscript of Pe denotes that it is obtained from either 
the temperature or velocity, and the typical velocity U is 
represented by j^;^ ■ As shown in Fig. [51 we can see that 
Pct and Pe„ have comparable magnitude, but their dif- 
ference can be well discerned. The Ra-dependence of Pex 
can be fitted with 0.55i?a"-^^. The exponent agrees excel- 
lently with the temperature oscillation measurements in 
low temperature helium gas 15 1, mercury [l^, and wa- 
ter [l^ of a cylindrical cell, the different prefactors are 
perhaps due to the fluids of different Pr. Moreover, the 
scaling is consistent with that of the angular oscillations 
of the horizontal thermal plume motion above the bottom 
plate in a cylindrical cell 2^. This result confirms that 



the temperature oscillation is the global and universal 
character in thermal convection, and this result also im- 
plies that temperature oscillation period (tx) is a result 
of flow motions of thermal structures and independent 
of system geometry. We also notice that some exper- 
iments have reported different scalings for temperature 
oscillation period (^t) in the cylindrical cell [l^ [j^]. 
One possible explanation is that the pathlength of the 
circulation affects the scaling of oscillation period [241] . 

In the cylindrical cell, it was reported that the veloc- 
ity oscillation follows 0.167i?a"'*^ and it is driven by the 
same mechanism as that of temperature 1^, 2^- That is, 
the velocity and temperature oscillations are generated 
by the alternating emission of thermal plumes between 
the upper and lower thermal boundary layers [lO] . More- 
over, it was found that both temperature and velocity os- 
cillation periods are coincident with the rotational period 
of the large scale circulation 17, 2^- In the present 



system, the i?a-dependence of Pe„ is fitted as 0.29Ra°-^^. 
In a previous PIV measurement [s^, the rotation rate 
of large scale circulation has been obtained, and the Pe 
could be fitted with the formula, 0.318i?a° ''^^, which is 
re-plotted in Fig. [Sj We can see the Peclet number Pe^ 
agrees this fitting line very well, but the temperature Per 
does not. These results suggest that the velocity oscil- 
lation period is consistent with the rotational period of 
the large scale flow in both cylindrical and rectangular 
cells, and both of them could be affected by the system 
geometry. Based on the results that temperature and ve- 
locity oscillations have different responses to the system 
geometry, it may be proposed that the two oscillations 
reflect different aspects of the thermal convective flow 
even though they are related to each other. 



CONCLUSION 

In this work, both temperature and velocity oscilla- 
tions can be found in the rectangular cell, where only one 
large-scale convection roll exists. At fixed Ra (8.9 x 10^^), 
it is found that temperature oscillation can be observed 
in most of the locations in which measurements were 
made, and they gave a position-independent oscillation 
period, = 74 ± 2 seconds. For measurements made 
with different Ra, the normalized oscillation period Pex 
is found to scale as 0.55Ra°'^'^ , which is the same as that 
found in cylindrical cells 0, El E i^. This scahng 
is also the same as that obtained from the horizontal 



motions of thermal structures (plumes) 2^ . This result 
suggests that temperature oscillation is a global and uni- 
versal characteristic of turbulent thermal convection, and 
its statistical behavior may be related to the flow motions 
of thermal structures and are independent of system ge- 
ometry. 

It is a surprise that the oscillation in the velocity field 
cannot be seen in its vertical component Vz near the side- 
wall. If the thermal bursts are erupted in a periodical 
way, it should be detected by the vertical velocity be- 
cause the thermal bursts are driven vertically under the 
buoyancy force. However, velocity oscillation can only 
be observed in its horizontal component (vy) at positions 
close to the sidewall, where its histogram exhibits two 
peaks. The normalized velocity oscillation period Pcy 
shows an Ra-dependence of Ra'^-^'^ and it collapses very 
well on that obtained from the rotation rate of convec- 
tive flow [3^ . This is different from the observation in the 
cylindrical cell where it was reported that the pe- 
riods for both the velocity and temperature oscillations 
and the rotational period of the convective flow follow 
the same scaling [13, [1^ . This comparison suggests that 
the velocity oscillation period is consistent with that ob- 
tained from rotation rate of the large scale flow, and both 
of them could be affected by the system geometry. 
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